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Starting from a quantum Langevin equation (QLE) of a charged particle coupled to a heat bath 
in the presence of an external magnetic field, we present a fully dynamical calculation of the sus- 
ceptibility tensor. We further evaluate the position autocorrelation function by using the Gibbs 
ensemble approach. This quantity is shown to be related to the imaginary part of the dynamical 
susceptibility, thereby validating the fluctuation-dissipation theorem in the context of dissipative 
diamagnetism. Finally we present an overview of coherence-to-decoherence transition in the realm 
of dissipative diamagnetism at zero temperature. The analysis underscores the importance of the 
details of the relevant physical quantity, as far as coherence to decoherence transition is concerned. 
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I. INTRODUCTION 



The problem of a quantum charged particle in the 
presence of a magnetic field is an old and important 
one 01 . When Landau gave the theory of diamagnetism 
a major breakthrough in solid state physics was made 
possible 0, Q- The physics of Landau levels is of 
great interest in the quantum Hall effect and high 
temperature superconductivity In the present 

paper we address the issue of what happens when we 
combine the Landau problem with the Drude transport 
treatment which naturally brings in the phenomenon of 
environment-induced dynamics |6|. 

The consequences of coupling of a system to its envi- 
ronment are threefold. First, energy may irreversibly be 
transferred from the system to the environment giving 
rise to dissipation [Si M H- Second, the spontaneous 
fluctuations injected by the environment into the system 
govern the response of the system degrees of freedom to 
weak, external stimuli [wL [ll| ■ Finally the entanglement 
between the system and environment degrees of freedom 
destroys the coherent superposition of quantum states, 
leading to decoherence [TJ • 

We discuss all the three above mentioned effects in the 
context of Landau diamagnetism which is inherently 
and intrinsically quantum in nature. For the purpose 
of investigating fluctuation, dissipation and decoherence 
in what we call Landau-Drude diamagnetism |l3j it is 
convenient to use the formulation given by Ford et al 

S UsI l , following the classical treatment due to Zwanzig 
Starting from the Feynmann- Vernon model in 
which a particle moving in an arbitrary potential is 
assumed to be linearly coupled with a collection of 
quantum harmonic oscillators these authors derive 
a quantum Langevin equation (QLE). We use this QLE 
as the basis of our further discussion, in what may 
be referred to as the Einstein approach to Statistical 
Physics [11. 

From QLE we calculate the dynamical susceptibility 
tensor. We then evaluate the position autocorrelation 
function from the Euclidean action of the Feynmann- 



Vernon model using the Gibbsian ensemble approach. 
This quantity is a measure of spontaneous fluctuations 
in the degrees of freedom of the system. We discuss the 
relation between the position autocorrelation function 
and the imaginary part of the susceptibility and thus 
establish the fluctuation-dissipation theorem in the 
context of dissipative diamagnetism. 
The destruction of quantum coherence by environment- 
induced dissipation is of central interest in atomic 
physics 0, condensed matter physics |23, as well 
as chemical and biological reactions We discuss 

this environment-induced decoherence in the context of 
dissipative diamagnetism. Landau diamagnetism has 
its origin in coherent circular motion of the electron in 
a plane normal to the magnetic field. This coherent 
motion is disturbed due to interaction with environ- 
mental degrees of freedom, e.g. defects, phonons, etc. 
We illustrate how the system transits from the coherent 
'Landau reg ime' to the decoherent 'Bohr- Van Leeuwen 
regime' [22, 123 ■ We show that the coherent-decoherent 
transition depends on the particular dynamical quantity 
(e.g., correlation function, occupation probability, etc.) 
under consideration. 

The paper is organized as follows. In the next section 
we discuss our model Hamiltonian and the correspond- 
ing QLE. In Sec. Ill we calculate the generalized 
susceptibility tensor. Sec. IV deals with the position 
autocorrelation function and its relation to the sus- 
ceptibility, thus establishing the fluctuation-dissipation 
theorem in the context of dissipative diamagnetism. In 
Sec. V we study the coherence-to-dccohercncc transi- 
tion. Finally we summarize our results and present a 
few concluding remarks in Sec. VI. 



II. MODEL, QLE AND EINSTEIN APPROACH 



We start with the Feynmann- Vernon Hamiltonian for 
a charged particle in a magnetic field B, coupled to an 



environment of quantum harmonic oscillators [1 
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where p and (f arc the momentum and position opera- 
tors of the particle, and A is the vector potential. The 
second term due to Darwin |24j represents a confining 
potential to recover the correct boundary contribution. 
Now following Ford et al IlilllSl one can write the QLE 
emanating from Eq.(l) as [13 



mci+ I dt'j{t~t')q{t') + mu;QCi~-{qxB) = F{t), (2) 

1 c 



where the auto-correlation and the commutator of F{t) 
are given by 

(^{Fo.{t),Fpit')}) = Sa,i3^ ^ 5R[7(a> + iO+)]huj 

cothi^^)cosMt-t')}duj, (3) 



([F^{t),Fp{t')]) ^6a0- / n[j{u; + iO+)]hu; 

\ ' Jo 

sm{Lo{t-t')}dLO, (4) 



where 7(5) = dtexp{ist)'-f{t) [Im s > 0). 
At this stage we introduce the nomenclature of Ohmic 
dissipation as well as non-Ohmic dissipation. Defining 
the spectral density of the environmental degrees of free- 



dom as J{uj) = § rrLjUj'jSiuj — ujj), we can rewrite 

the memory kernel "f{t) in terms of the spectral density 
as 
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where Q{t) is the Heaviside step function. In the 
strictly Ohmic case, damping is frequency-independent 
and the spectral density J(ti^) = rwyuj. The memory 
kernel ^{t — t') is thus replaced by mjS{t — t'), so that 
3?[7(a; + iO^)] reduces to 7717, a constant. In this limit 
we get an ordinary Langevin equation. It is interesting 
to note that the underlying stochastic process is still 
non-Mar kovian, even though there is no memory. In the 
non-Ohmic case (e.g. due to a phonon bath) the spectral 
density is defined as J(w) = 7717(0;), where 7(0;) = 70;^. 
The damping kernel 7(1^) then brings in memory-friction 
effects. 

III. GENERALIZED SUSCEPTIBILITY 
TENSOR 

In this section we consider the linear response of the 
position coordinate to an external force /(i), assumed 
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FIG. 1: (color online) The imaginary part of susceptibility 
Xxx (a) Ohmic dissipation case ( J(a;) ~ lj) for two ujc values, 
(b) Ohmic dissipation case for two 7 values, (c) Non-Ohmic 
dissipation case {J{(^) ~ OJ^) for two uic values, (d) Non- 
Ohmic dissipation case for two 7 values. 



small. By imagining the force to have been switched on at 
time t = — 00 all transient effects can be ignored and the 
nontransient response can be captured by the frequency- 
dependent generalised susceptibility. The corresponding 
QLE now reads 

•• /■* • e • ^ 

mq+ / dt'-f{t - t')q{t') + mujlq ((fx B) = 

J-oa c 

Fit) + fit). (6) 
We rewrite Eq. (6) in a Fourier transformed form 

imiiol - cj2) - iuj^iuj))Sai3 + iuj-eafjpBp qpiuj) = 

Fo.iLu) + UiLu), (7) 

with 
Ziiuj) = 



dte^'^^Z.it) 

(i = 1,2,3,4; Zi = 7, Z2 = qp, Z3 ^ F^, Z4 = /„), 



and eat3p is the Levi-Civita symbol, a, /?, p being the three 
spatial directions (i.e. a, /3, p = x, y, z). In linear response 
theory one can write |25| 



Qait) 



dsXc.pit-s)iFpis)+ffiis)), (8) 



where Xq/3 is the generalised susceptibility tensor. In 
Fourier transformed form Eq. (8) becomes 

9a(w) = Xa/3(^)[F/3(^) + M^)]- (9) 

Rewriting Eq. (7) as 

Yo,i3icj)q0iu;) = [Fo^ioj) + Uiu;)], (10) 
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FIG. 2: (color online) The real part of susceptibility Xxx (a) 
Ohmic dissipation case (J(tt') ~ uj) for two Uc values, (b) 
Ohmic dissipation case for two 7 values, (c) Non-Ohmic dis- 
sipation case {J{co) ~ w^) for two lOc values, (d) Non-Ohmic 
dissipation case for two 7 values. 



the generahsed susceptibiUty can be evaluated as 

X0./3 = [Y-\u;)U, (11) 
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where A{ll>) — m(ti>g — lu'^) — iuj^{uj). Clearly 
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where 
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where (*) denotes the complex conjugate of the corre- 
sponding variable. 

The expression is simplified when the magnetic field is 
taken along z axis, thus 
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For this particular case the real part of the susceptibility 
is 
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where the cyclotron frequency ojr = — . For the Ohmic 
dissipation case the susceptibility has four poles at 



(jj = ±UJ+ = [ ± — — J 



^ = = ['"-+^^ ± vM±^EZE2^^](i8) 

On the other hand, for the non-Ohmic case these 
poles cannot be evaluated analytically. The numerical 
results for the Ohmic dissipation as well as non-Ohmic 
dissipation cases are presented below. 
We plot in Fig. (1) the dissipative part of the x- 
componcnt of susceptibility i.e. Xxxi^) versus oj for 
different values of cUc and 7 in accordance with Eq. 
(17). We note that Xxxi^) is odd in lu for the Ohmic 
dissipation case and has Lorentzian line shapes for finite 
damping values, with peaks centered at the poles. For 
the non-Ohmic case Xxxi^) even in uj. It is evident 



4 



from Fig. (1-b) that for finite but weak damping one can 
obtain all the four peaks for the ohmic dissipation case 
whereas for high damping only two peaks are obtained. 
The same is true for the non-Ohmic case (Fig. (1-d)). 
The only difference is that the magnitude of the peak 
height is higher for the non-ohmic case and is always 
positive. Also the peak width increases with the increase 
of 7 for both Ohmic and non-Ohmic cases. On the 
other hand the width of the peak decreases with the 
increase of Wc, as is expected on physical grounds. In the 
non-Ohmic case the number of peaks also increases from 
two to four with the increase of Wc, whereas it remains 
two for the Ohmic case with the increase of if 7 is 
kept large. Thus, dissipative effects are stronger for the 
Ohmic case. 

In Fig. (2) we plot the reactive part or the real part 
of the x-componcnt of susceptibility {x'xxi'-^)) versus 
u for different values of ujc and 7 in accordance with 
Eq. (16). x'xxi^) is odd in to for the Ohmic as well as 
non-Ohmic cases. The spreading of the peaks increases 
but the peak height decreases with the decrease of Uc for 
the Ohmic case. On the other hand both the spreading 
and peak height decrease with the decrease of uJc for the 
non-Ohmic case. But the features are same with the 
variation of 7 for both Ohmic and non-Ohmic cases — - 
the peak height increases but the spreading decreases 
with the decrease of 7. In addition the number of peaks 
increases from one to two with the decrease of 7 in the 
Ohmic as well as non-Ohmic cases. 

The z component of the susceptibility tensor is of course 
the same as that of a damped harmonic oscillator 
because it has nothing to do with ujc- 

IV. FLUCTUATION - DISSIPATION 
RELATIONSHIP: GIBBS APPROACH 

In Sec. Ill we calculated the susceptibility as the 
asymptotic (i.e. t — > 00) response from a fully time- 
dependent formulation of the underlying QLE. Because 
detailed balance relations (viz. Eqs. (3) and (4)) are 
built-in within the QLE, as the heat bath is assumed to 
be in thermal equilibrium at a fixed temperature T, the 
asymptotic response is expected to be related to the equi- 
librium properties of the system. This expectation is at 
the heart of what Kadanoff calls the Einstein approach to 
Statistical Mechanics |l8l | in which equilibrium answers 
are sought to be obtained from the asymptotic limit of 
time-dependent results. It is then natural to ask whether 
the response obtained from the Einstein approach can be 
related to spontaneous or equilibrium fluctuations that 
can be independently calculated from the standard Gibb- 
sian formulation of equilibrium Statistical Mechanics. If 
we can establish this relation it will not only be tan- 
tamount to establishing the fluctuation-dissipation the- 
orem for the phenomena at hand, but also to demon- 
strating the equivalence of the Einstein and the Gibbs 
approaches to Statistical Mechanics |2^ . 
With this preamble the position autocorrelation function 



in equilibrium is defined as 

C{t) =< x{t) ■ f(0) >= Tr{x{t) ■ x{0)pi3), (19) 

where pp is the equilibrium density matrix of the full 
system and x is the two dimensional position vector in 
the x-y plane.. We determine C{t) by first calculating 
its imaginary time version starting from the Euclidean 
action of the system as described by Eq. (1) 
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where the first term takes care of the system part, the 
second term accounts for the coupling to the environ- 
ment and the third term corresponds to the interaction 
with an external force in imaginary time. This helps us 
to determine the correlation function by variation with 
respect to this force [27ll28| 



< x(t) • x{a) >= h^Tr 
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It is sufficient to restrict ourselves to the classical 
path for the calculation of the autocorrelation function 
[271 123 ■ Thus the Fourier representation of the classical 
Euclidean action becomes Ha, 123 
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where I'n = are the so-called Matsubara frequencies. 
Since the force appears only through the action in the 
exponent of the equilibrium density matrix, we can eas- 
ily evaluate the functional derivatives according to Eq. 





FIG. 3: The analytic continuation of the imaginary time cor- 
relation function to real times by using the contours depicted 
in (a) and (b) to obtain Eq. (24) and Eq. (25) respectively. 
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(21) and obtain the position autocorrelation function in 
imaginary time: 



1 



exp(w„r). (23) 



by simply replacing r by it, because for negative times 
the sum does not converge. The idea is to express the 
sum in Eq. (23) as a contour integral in the complex 
frequency plane. To do this, we need a standard trick of 
Statistical Mechanics 12^. Look for a function which is 
well-behaved at infinity, but has poles at uj = ivn- This 
requirement is fulfilled by . Now doing the 

integration along the contour shown in Fig. (3-a) we find 



The real time correlation function cannot be obtained 
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Similarly, an integration along the contour shown in Fig. (3-b) leads to 
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Subtracting Eq. (25) from Eq. (24) we obtain 
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We may now pass to real time by the replacement t it to obtain the real time correlation function 
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It is easy to show from Eq. (27) that 



1 — exp[—(}fko) 

Eq. (28) represents the fluctuation-dissipation theorem 
in the context of dissipative Landau diamagnctism. The 
position autocorrelation function describes the sponta- 
neous fluctuations of the system while the imaginary 
part of the dynamic susceptibility Xxx determines the 
energy dissipation in the system due to work done by an 
external weak force. 



V. COHERENCE - DECOHERENCE 
TRANSITION 

In this Section our discussion is focused on the destruc- 
tion of quantum coherence by environment-induced dis- 
sipation in the context of dissipative diamagnctism. Two 
questions arc relevant: (i) Can we quantify the criterion 
for crossover from coherent to decoherent dynamics? (ii) 
Is this criterion universal? As far as some model systems 
are concerned, the answer to (i) is in the affirmative 33 ■ 
Regarding the question (ii), there seems to be no uni- 
versality in the criterion of crossover. As a matter of 
fact, the value of the crossover parameter depends on 
the particular quantity under consideration and its ini- 
tial preparation. Thus, quantum memory effects play a 



crucial role as the system makes a transition from the 
coherent to the decoherent regime. To clarify this issue 
we focus on dissipative diamagnetism and we consider its 
T = behavior, wherein quantum coherence is the most 
prominent. Here we follow the discussion of Egger et al 

We start with the QLE for dissipative Landau diamag- 
netism subject to Ohmic damping. The motion in the 
x-y plane can be expressed in the compact form: 



'7 I -7 I 2^ ^(*) 
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where Z = x + iy, j = j + iujc, and 9 = 
Thus, the time-dependence of the corresponding classical 
quantity {a la Ehrcnfcst) is governed by the following 
equation: 



< Z > +7 < Z > +L0I < Z >= 



m ' 



(30) 



where the angular brackets represent statistical averages 
over the ground state properties (T = 0), i.e. the ex- 
pectation values. As discussed earlier the response to an 
external force is characterized by the generalized suscep- 
tibility Xoscit) M- 



< Z{t) >= 
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dt' Xoscit - t')e{t'). 



(31) 



From Eqs. (24) and (25), we obtain the Fourier trans- 
form of Xosc{t) as 



Wo 



(32) 



On the other hand, using the fluctuation-dissipation the- 
orem Xosc(w) can be related to the spectral function 
Sosc{^), which in turn determines the equilibrium cor- 
relation function Cosc{oj)- The functional relationship 
which holds at T = is as follows: 
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FIG. 4: (color online) Spectral function Sosc{lo) vs. ui with 
Ohmic dissipation for dissipative Landau diamagnetism for 
different parameter values. 



are evident from Fig. (4) in which we plot Sosc{^) vs. 
Lj for different ki and ^2- At the critical coherent cri- 
terium (c.f. Eq. (37)) the two peaks merge into a single 
quasielastic peak. The latter persists for > k^. Since 
the quasielastic peak is centered near u; ~ 0, we can make 
a small-o; expansion of Sosc{^)'- 



Soscioj) ^ KlXo 1 ~ 1^2X0^ + (2 - - K2)Xo'^^ + 0{UJ^) 

(35) 

where Xo = — • The critical line is determined by in- 
specting the sign of the curvature of Sosci^)- The latter 
is positive (implying coherence) if ^^^'2''^^ > 0, or 



\uj\ 

where Cosc{t) = 3? < Z{t)Z{Q) >. Using Eqs. (32) and 
(33), we obtain the spectral function 



(34) 



The quantity Soscii^) can be used as a signature for the 
transition from coherence to decoherence: Soscii^) has 

two inelastic peaks at ~ ^ — K2 ± ■\/4 — + 
for weak damping, where ki and K2 are dimcnsionlcss 
parameters defined by ki = and K2 = These two 
quantities are employed as the crossover parameters to 
quantify the coherence to decoherence transition. Defin- 
ing = k1 + K2, we can say that below the critical coher- 
ent criterium (defined below, c.L Eq. (37)) i.e. < k^, 
the function Soscii-^) exhibits two inelastic peaks which 



/s: — 2. 



(36) 



But the curvature changes sign at the critical line: 

rI^kI + kI^ 2, (37) 
and hence the system goes to the decoherent region when 



Kn > 2. 



(38) 



It is illustrative to compare this behavior with that of 
the damped harmonic oscillator which was discussed by 
Egger et al HO]. 

Comparing Figs. (4) and (5) one notes that for the 
damped oscillator case Sosc{^) has two inelastic peaks of 
equal height for weak damping. As the one-parameter 
damping strength increases these two peaks approach 
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peak persists. 

We turn next to a different criterion for quantifying the 
transition from coherence to decoherence, which is based 
on the quantity Posc{t), defined as follows: 



< Z{t) > 



(39) 



We are interested in the relaxation of the expectation 
value < Z(t) > starting from a noncquilibrium initial 
state. Applying the force F{t) = mt^pZo for t < 0, the 
initial condition < Z{0) >— Zq is prepared and the cor- 
responding dynamical quantity Posc{t) is computed, af- 
ter switching off the force F{t), at t = 0. Following the 
damped quantum harmonic oscillator case [3(1 | we may 
now write 



Poscit) = 



where 



cos{flt — (j)) exp(— ^) 
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FIG. 5: (color online) Spectral function Sosc{(j->) vs. u with 
Ohmic dissipation for a damped harmonic oscilltor for differ- 
ent parameter values. 



each other and at the critical damping strength (ac) the 
two peaks merge into a single quasielastic peak at w = 
which persists for a > ac- On the other hand, for dissi- 
pative diamagnetism, the coherent-decoherent transition 
is to be examined in a two-parameter plane, defined by 
Ki and K2- One obtains two inelastic peaks which are not 
of equal height for low values of ni and K2 because the 
peaks are not symmetric in either side of w = 0. As one 
increases ki and K2 the peak height of the small peak 
decreases and at the critical line the small peak vanishes 
and we obtain a single peak which is not at w = 0, but 
near tj = 0. Above the critical line the single quasielastic 
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The signature of coherence is now damped-oscillatory be- 
havior if 6^ > and + iP' > 0. Thus the important 
inequality condition for the system to be coherent is: 



(1 



(44) 



The system crosses over to relaxational (decoherent) be- 
havior at the critical line 



= 0, 



(45) 



which is clearly different from the criterium mentioned 
above (cf. Eq. (37)). Thus the criterion for crossover 



from coherence to decoherence depends on the specific 
physical quantity considered. This conclusion is similar 
to the cases of damped quantum harmonic oscillator as 
well as the spin-Boson model ^i^. 



VI. SUMMARY AND CONCLUSION 

We have analyzed here an exact treatment of the 
Feynmann- Vernon model of a charged Brownian particle 
in a magnetic field in the quantum dissipativc regime. 
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Starting from the QLE wc have derived the generahsed 
susceptibility tensor, and have discussed its real and 
imaginary parts for the particular case when the mag- 
netic field B is along the z axis. Following the Gibbs 
ensemble approach, we have calculated the position 
autocorrelation function that measures the spontaneous 
fluctuations of the system degrees of freedom due to 
coupling with the environment. The latter has been 
shown to be related to the imaginary part of the sus- 
ceptibility that measures the energy dissipation of the 
system due to irreversible energy transfer between the 
system and the environment. The aforesaid treatment 
then exemplifies the fluctuation-dissipation theorem 
in the context of dissipative diamagnetism as well as 
proves the equivalence of the Einstein and the Gibbs ap- 
proaches to Statistical Mechanics. Environment-induced 
decoherence is an important issue in mesoscopic systems 
and quantum information processes. We have discussed 
this in the context of dissipative diamagnetism and 



have argued that the transition from the Landau to the 
Bohr- Van Leeuwen regime can indeed be viewed as a 
coherence to decoherence transition. Further it has been 
demonstrated that the initial preparation of a dissipative 
quantum system leads to abrupt changes regarding the 
criterion for coherent to decoherent transition. As in 
glassy systems characterized by hysteretic behavior, 
quantum systems also exhibit memory of their initial 
state of preparation. In conclusion, we have presented a 
unified treatment of threefold response, i.e. fluctuation, 
dissipation and decoherence of a system due to its 
coupling with environment in the context of the con- 
temporarily important topic of dissipative diamagnetism. 
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